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I ' The heart rate variability (HRV) in diabetic human subjects, has been analyzed using 

^ '. 
1^ , lagged Poincare plot, auto-correlation and the detrended fluctuation analysis methods. 

The parameters SDl, and SD12(= SD1/SD2) for Poincare plot for diabetic are lower 

than that for non-diabetic subjects and reverse is case for SD2 for all lagged number (m). 

c/3 ■ The slope and the curvature of the plot SD12 vs m is much reduced for diabetic subject. 

jj^ ' The scatter plot of two successive interbeat intervals points out smaller variability in 

o ,' 

diabetic heart. The detrended fluctuation exponent has a higher value for diabetic group. 

The auto-correlation function of the deviation of interbeat interval in diabetic group shows 

highly correlated pattern when compared to that of normal one. The study suggests that 

^^ ' the curvature of SD12 and auto-correlation method appear to be better way to assess the 

lO I alteration of regulatory system on heart dynamics in diabetic condition. 
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1 Introduction 

The variability of human heart rate is well documented and unraveling of its complexity has 
drawn continuing attention. The cardiac oscillators are embedded in dynamic environment and 
are under the influence of regulatory biological system like autonomic nervous system (ANS). In 
order to cope with influencing signals some amount of plasticity in cardiac activity are therefore 
essential and this generates the variability of heart rate intervals. Heart Rate variability (HRV) 
is a measure of the fluctuations in the intervals between heartbeats, the RR interval. It reflects 
the cardiac autonomic regulation [1-3]. Diabetes, a common disease, may lead to the serious 
complication of autonomic diabetic neuropathy. Autonomic diabetic neuropathy is a major risk 
factor for cardiovascular mortality [4] . It leads to impaired regulation of blood pressure, heart 
rate and heart rate variability (HRV). It is associated with increased rate of cardiovascular 
disease such as heart attack, sudden cardiac death or Silent ischemia [5-9]. Early subclinical 
diagnosis of Autonomic Diabetic Neuropathy is difficult and may not be very sensitive. One 
of the common methods of detection is the Ewing Test Battery. Ewing Test Battery consists 
of flve tests and requires high cooperation from the patient. Its use is also limited by the fact 
that it cannot be used in patients with co-morbidities like existing heart disease or respiratory 
disease, which are contraindication for the Valsalva manoeuvre, an integrated part of the battery 
[10]. Newer non invasive approaches like spectral analysis and Heart rate variability analysis are 
completely non invasive and are not associated with such limitations. They hold promise to 
provide a better diagnostic or prognostic tool for Diabetic Autonomic Neuropathy [11-14]. 

However, one has to be cautious while using the heart rate variability as conventionally used 
time and frequency domain parameters of HRV may not always represent the non- stationary 
characteristic of EGG rooted in nonlinear dynamics of the heart. A non linear method should 
be an appropriate agent to analyze the heart rate variability. Poincare or Lorenz plot, a scatter 
plot of each RR interval as a function of the preceding one, is a non linear method with the 
ability to capture the non linear pattern of heart rate dynamics [15-18]. A recent approach 
with auto correlation technique also provides a new tool for this non linear system analysis 
[19]. In this study the time series of RRn -interval for number of diabetic and control subjects 
have been analyzed following the Poincare plot and its modifications. Aim has been to reveal 
relative importance of the plot-characteristic parameters in differentiating a diabetic heart from 
non-diabetic one. We also stress that this investigation is of exploratory nature to determine 
the infiuence of the regulatory system on diabetic heart as elucidated by different kind of 
analysis of RR-interval. A marked difference between a diabetic and normal heart dynamics 
has been observed in SD12, correlation between successive difference in RRn -interval and the 
auto- correlation of fluctuation of RR„ . 



2 Method 

We selected 23 patients to carry out our research. We applied a standard exclusion criterion 
to ensure that any changes in HRV detected were due to the effect of the diabetes. A group of 
23 subjects was assigned as the control after matching the age distribution and other probable 
confounding parameters. We have maintained the principles outlined in the Declaration of 
Helsinki and a proper informed consent was given by every subject. All subjects were instructed 
to avoid caffeine, alcohol and heavy exercise the day before the study. The subjects were asked 
to relax for a period of 15 minutes before we took the recording for 10 minutes in supine 
position. One limitation of the study population was the non uniformity of medication among 
the study population. 

3 Electrocardiogram (ECG) recording and HRV measurement 

Three limb electrodes were placed on the appropriate limb to obtain the ECG signal. The 
sampling rate was lOOOHz. Using A/D converter the ECG signal was recorded into a computer 
for analysis. Ectopic beats (if any) were selected visually and deleted manually from recorded 
data for 10 minute interval. Then, RR intervals (in fact, NN intervals, i.e., normal-to-normal 
intervals) were measured and acquired through software written in Origin and Matlab. 

4 Poincare plot indices: 

Poincare plot, as described earlier, is basically a scatter plot of RRn and RRn+i- When the 
plot is adjusted with an ellipse, the analysis provides three indexes: S'-Dl, SD2 and ratio 
S'-D12 = SD1/SD2 [16]. SDl is the standard deviation of instantaneous beat to beat interval 
variability and it is width of the ellipse and SD2 is the length of the ellipse. It has been 
shown that S'-Dl correlates with the short term variability of the heart rate and SD2 is a 
measurement of long term variability [17]. A strong correlation has been established between 
the high frequency powers of Heart rate signals (modulated by parasympathetic nervous system) 
to SD2. A similar correlation has been found between SD2 and both low and high frequency 
power (modulated by both the parasympathetic and sympathetic nervous system) [20]. In 
addition to this conventional plot {RRn+i vs RRn) we also used the generalized Poincare plot 
with different intervals i.e the m lagged Poincare plot ( plot of RRn+m against RRn)- The 
concept of this m lagged plot emerged from the recognized notion that any given RR interval 
can influence up to eight subsequent RR intervals [21, 22]. 



5 Autocorrelation technique 

Application of autocorrelation technique is a very recent idea and has enough potential to 
provide a new insight towards autonomic regulation of heart [19]. In this approach HRV is 
considered as the interaction between coupled oscillators of various frequencies [23]. It also 
reflects on the imbedded time scale of HRV. It has been taken as each of these time scale in the 
coupled oscillator is represented by a separate self-oscillator interacting with other oscillators 
with a particular physiological function. It enables us to count for the contribution of various 
systems like respiration, neurogenic, myogenic and endothelial towards the HRV. In this study 
we have applied the correlation technique to RR interval through the Origin software. 

6 Results and Discussions 

The plot of RRn+m with RRn for m = 1, 5 and 9 is depicted in Fig.l for two subjects of similar 
age but one is diabetic (left figure) and other is non diabetic and disease free (right figure). 
With an increase in the lag number the plot becomes more scattered with consequent incre- 
ment in both width and length of the plot. For diabetic subject the width of the plot RRn+m 
vs RRn is smaller compared to that for non-diabetic person. SDl and SD2 were calculated 
for lag m from the relation :SD1 = (<l>(m) - $(0))^/^ and SD2 = (<l>(m) + $(0))^/^ where the 
auto-covariance function <l>(m) is given by 

$(M) =E[{RRn - R:R){RRn+m - RR)]- 

where RR is the mean i?i?„ [17]. The extent of scatter in lagged Poincare plot is quantified with 
the measure SDl, SD2 (both in sec.) and SD12, and the variation of these measure with lag 
number m is presented in Fig. 2 (points). Each point in the plot 

is obtained from calculation of the above parameters for individual subjects and then by 
taking the mean of that group. The values of SDl and SD12 are higher in non-diabetic 
group compared to diabetic one, and the situation is reversed in case of SD2. The difference 
between the values of SD12 (points) in normal and diabetic subjects are more significant for 
different lag number. The differences between the values of all these parameters for two groups 
are statistically significant with a p value less than 0.001. In order to find the analytical 
relationship between S'-D's and m we used the method of Fade Approximant [24]. We assumed 
a simple form of the Fade approximant for SD^s as the ratio of polynomial in m of degree one. 

a + bm 1 + I3m .. 

^ = I J = x-r, — (1) 

c + am 1 + '-)m 



Table 1: The value of parameters X) P) 7 obtained from fit of eq.(l) with respective value of 
R^ .The parameters L and Q are the coefficient of linear and quadratic terms in expansion of Y 
in terms of m. Values of x, L and Q for SDl and SD2 are in second. 





a 


XxlO-2 


/3xlO-2 


7x10-2 


i?2 xlO-2 


L xlO-3 


Q xlO-4 


SDl 


ND 


1.3 ±0.03 


39.1 ±2.0 


3.2 ±0.2 


99.9 


4.7 ±0.4 


1.5 ±0.2 




D 


1.00 ±0.02 


38.2 ±1.4 


2.00 ±0.1 


99.9 


3.6 ±0.08 


0.7 ±0.02 


SD2 


ND 


3.2 ±0.06 


20.3 ±1.1 


3.5 ±0.2 


99.9 


5.4 ± 0.04 


1.9 ±0.2 




D 


3.1 ±0.07 


26.4 ±1.6 


4.4 ±0.3 


99.9 


6.8 ±0.6 


3.0 ±0.5 


SD12 


ND 


40.2 ±0.5 


25.0 ±1.8 


12.2 ±0.9 


99.9 


51.3 ±6.4 


62.7 ±2.4 




D 


33 ±0.3 


15.3 ±.8 


6.5 ±0.4 


99.9 


29 ±2.6 


18.9 ±2.8 



Here Y = SD1,SD2 or SD12 and x = '^/c, /? = b/a and 7 = d/c were taken as new 
unknown parameters. An excellent fitting of the data with the equation (1) (solid line on the 
curve) was found with different set of x, P, 7 which were listed in table-1, and in all cases R^ 
value was 0.999. In order to find relevance of these parameters in assessing health of cardio- 
vascular system we considered the eq. (1) for small m. In this limit the equation (1) can be 
approximated as Y = C + Lm ± Qm^ where the slope L = x(/3 — 7) and the curvature Q = - 
7L. So the slope and curvature of plot of SD vs m were determined by the fitted parameters x, 
P, 7. Earlier, it was noted that the curvature of the plot {SD12-m) was significantly decreased 
for subject with cardio- vascular disease [21]. The values for L and Q as obtained from the fit 
of data by eq. (1) were given in table-1. The general features are that the slope is positive but 
curvature is negative for all parameters and curvature is nearly one order of magnitude smaller 
than the slope. 

The cardiac oscillator is embedded in the dynamic environment regulated by ANS. This 
in turn leads to the correlation between interbeat variability. The coherent nature of RRn 
interval can be assessed from the map of interval variation: rr„ . 1 = — "^pT? ^""""^ vs rrn = 
— "j^^- " where \RRni represents the mean of RRn- This plot reveals the correlation between 
the variabilities of three consecutive RR intervals. Without any modulating effect the plot 
reduces to a point at origin whereas in presence of random input uniform distribution of point 
is expected. In Fig .3 the map of rr^+i -rVn for two subjects from each group (ND -control and 
D - diabetic) has been shown. The subjects of two groups are age matched. 

Distinct differences are evident in the map. For diabetic patients the points are more 
crowded around origin. For diabetic person (Fig. 3 Dl) two successive decelerations (first quad- 
rant) and acceleration (third quadrant) are as probable as combination of one acceleration and 
one deceleration (second and fourth quadrants). For D2 small asymmetry with more points 



in second and fourth quadrants is observed. In contrast, there exits more scatter about origin 
and higher asymmetry for normal person. Considering similar analysis for more subjects, the 
general features that emerge are high density of points around origin and smaller asymmetry 
for diabetic and more scattered with larger values of rVn and higher asymmetry for control 
group. It suggests that the influence of decelerating and accelerating mechanisms that act in 
different time frame on cardio-dynamics is substantially reduced in diabetic patients. 

The reduced value of short term beat fluctuation SDl indicates that the parasympathetic 
regulation in diabetic heart is weakened as expected. The higher value of SD2 in diabetic 
group means increased long term variability. Combination of these leads to higher ratio of 
SD12 in control group. The results for the Poincare plot are more pointed when one examines 
the values of the slope (L) and curvature (- Q) of the plot. For diabetic group, L and -Q for 
SDl and SD12 are smaller but their values for SD2 are higher compared to corresponding 
results for control group. The change in Q between two groups is more significant than the 
corresponding change in the values of L. In particular, the value of the curvature (Q) of SD12 
in control group is found to be more than three times of Q for diabetic group. It may be 
mentioned that very small value of curvature is found in patient with cardio-vascular disease 
(21). Normally, SD2 is well correlated with both high (HF) and low (LF) frequency components 
of RRn fluctuation spectrum and SDl is associated with HF.This data strongly suggest the 
reduction of parasympathetic activity and the over influence of sympathetic activity in the 
diabetic heart. Interestingly it also indirectly proves that sympathetic over influence has a 
strong correlation with cardiac morbidity [26, 27]. An increased 5*1^12 is considered as good 
indicator of heart dynamics. The study supports notion of reduced sympatho-vagal balance in 
diabetes. 

The strength of the correlation of heart rhythm can also be estimated by considering the 
auto-correlation of the deviation oi RRn from its mean < RRn >■ The auto-correlation function 
C{m) of a particular subject was calculated from C{m) = "^n — "^"^ — - ■ Where the 
deviation ARRn = RRn— < RRn > and N is total number of RRn interval. For representative 
result same two persons from two groups were considered and the results are displayed in Fig. 4. 
The results for two diabetic subjects (Dl and D2 of Fig 3) and for two normal subjects (NDl 
and ND2 of Fig.3) have been presented in upper and lower figures respectively. There also exits 
a distinctive feature in auto-correlation function for these two groups. For diabetic 

subjects, the correlation function C{m) decreases slowly (black and green curve of upper 
figure) with lag time. The time dependence was close to the sum of the two exponentials with 
superimposed small amplitude oscillation of low frequency. The presence and importance of 
the correlated pattern of beat interval was assessed by randomizing the RR-interval of each 
subject. The actual time series of RR- interval was shuffled using Matlab programme and 



the auto- correlation function of shuffled data (red and blue) were plotted in same Fig. 4. The 
function has very small correlation time that characterizes the exponential decay of C{m) . 
This indicates that a long time memory effect exists in diabetic condition. On the other hand 
the function C{m) for actual time series of normal persons are found to have exponential fall 
with smaller correlation time compared to that for diabetic cases. The correlation function for 
normal subject could be fitted with form C{m) ^ exp{—m/r)[a + bcos{ujm)]. However, the 
auto-correlation functions of shuffled data of all subjects (2-diabetics and 2-non-diabetics) were 
nearly identical. The memory time ( Tx < RR >) was much smaller compared to the diabetic 
subjects. The properties of ARRn can also be characterized by the probability distribution 
function P (ARRn). In (Fig.5) the result of P (ARRn) is given for Dl and 

NDl subject. For diabetic person the probability distribution is almost symmetrical and 
fits with a Gaussian fit (-R^ = 0.93) (full curve) with width = 0.023. For normal person the 
probability distribution P is asymmetrical with positive mean and higher width = 0.036 as 
obtained by the Gaussian fit (-R^ = 0.93). 

Another way to assess long term correlation in RR-time sequence is based on the Detrended 
Fluctuation Analysis (DFA) [25]. The measure of correlation was given by a scaling exponent 
(a) of the fiuctuation function F{t) ^ r" . The computation of fiuctuation function F{t) was 
done in the following way. For a given time sequence R{ti), ti = i6t where 6t is characteristic 
time interval for the sequence and i=l,N , an integrated time series r(tj) was defined as r(tj) = 
J2)[R{'tj)~ < RR >]) "^ = 1?^ • The integrated series was divided into boxes of equal size of 
time r = n6t and linear function was used to fit box data. The fiuctuation function F{r) was 
calculated as root mean square fiuctuations relative to the linear trend and a was obtained 
from fitting power law behaviour. It has been observed that acceptable estimate of the scaling 
exponent a (from DFA) can be obtained from analysis of data sets of length 256 samples or 
greater (equivalent to approximately 3.5 min for RR data at a heart rate of 70 bpm). The 
analysis of RR data for period of 10 min time interval was therefore expected to provide an 
adequate measure of the scaling exponent. The results are depicted in Fig. 6 for two groups. 
The mean value of a = 0.88pm. 17 for normal person is smaller than that a = 1.02 ± .13 for 
diabetic with p < 0.001. 

7 Conclusions 

Diabetic condition is marked by much lower value of SD12 for all lag number. The data 
on parameters SDl, SD2 and SD12 can be mapped by simple formula. The parameters 
have same functional dependence on lag number m suggesting similar source of variability of 
these quantities. The drastic reduction of the curvature of SD12 points out its significance in 



determining influence of regulatory system and environment and thus the health of the heart. 
The lagged Poincare plot of RR-interval appears to be better suited for assessing heart rate 
variability. Among the other analysis methods we have applied, we would like to emphasize 
on the novel application of autocorrelation technique to assess the autonomic regulation of a 
diabetic heart. To our knowledge this is the first attempt to differentiate between a normal and 
diabetic heart using autocorrelation technique. Results from the autocorrelation analysis of 
RRn interval suggest lack of short term variability of RR intervals among the diabetic patients. 
On the other hand, higher short term variability in normal subjects is marked by the presence of 
appropriate oscillation in correlation function. Diabetes as hypothesized impairs this variation 
and changes the RR interval distribution pattern. Moreover these change can be pictorially 
expressed in such a manner that even an untrained eye can distinguish the marked difference. 
These may prove to be really effective in an early diagnosis of diabetic autonomic neuropathy 
and may also have some prognostic importance. 
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Figure 1: Poincare plot RRn+m vs -R-R„ of one diabetic(D) and one non-diabetic (ND) subject 
of similar age for M = 1, 5 and 9. 
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Figure 2: Variation of mean 5/^1 (upper), mean SD2 (middle) and mean SD12 (lower) with lag 
number m for diabetic (D) and non-diabetic (ND) groups 



0.15^ 



o.m 



0.05^ 



S o.oa 



-0.05 



-0.10 



-0.15 




-0.15 -0.10 -0.05 0.00 0.05 0.10 0.15 

rr 

n 



0.15 



0.10 



0.05 



= 0.00 



-0.05 



-0.10 



-0.15 



D1 






^ 1 ' 1 ' 1 

■0.15 -0.10 -0.05 0.00 0.05 0.10 0.15 



rr 





0.15, 


I ND2 


I 


I 




0.10 
0.05 


1 


■ ■ 

.■■■■ 




'T 


oo 


Ji 








M 






-0.05 


^ffl 


7' 






-0.10 


1 ■ 

I 


■ 






-0.15 


1 1 1 1 1 


' 1 ' 


' 



0.15 



-0.15 -0.10 -0.05 0.00 0.05 0.10 0.15 

rr 



0.10 



0.05 



0.00 



-0.05 



-0.1 



-0.15 



D2 




^ I ' I ' I 

0.15 -0.10 -0.05 0.00 0.05 0.10 0.15 



rr 



Figure 3: Plot rr„+i vs rvn of diabetic(D) and non-diabetic (ND) subjects of similar age. 




Figure 4: Plot C{m) vs m of two diabetic(Dl and D2 ) and two non-diabetic (NDl and ND2) 
subjects considered in Fig. 3 
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Figure 5: Plot of probability P as a function of rvn of diabetic(Dl) and non-diabetic (NDl) 
subject of Fig. 3 
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Figure 6: Exponent a for non-diabetic (ND) and diabetic(D) groups 



